Na 2 Ir03 as a Novel Relativistic Mott Insulator with a 340 meV Gap 
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We study Na2lr03 by ARPES, optics, and band structure calculations in the local-density approx- 
imation (LDA). The weak dispersion of the Ir 5d-t2 g manifold highlights the importance of structural 
distortions and spin-orbit coupling (SO) in driving the system closer to a Mott transition. We detect 
an insulating gap A gap — 340 meV which, at variance with a Slater-type description, is already open 
at 300 K and does not show significant temperature dependence even across Tiv — 15 K. An LDA 
analysis with the inclusion of SO and Coulomb repulsion U reveals that, while the prodromes of an 
underlying insulating state are already found in LDA+SO, the correct gap magnitude can only be 
reproduced by LDA+SO+U, with (7 = 3eV. This establishes Na2lr03 as a novel type of Mott-like 
correlated insulator in which Coulomb and relativistic effects have to be treated on an equal footing. 

PACS numbers: 71.20.Be, 74.25. Jb, 74.25.Gz, 71.15.Mb 



The proposal of an effective J e // = l/2 Mott- Hubbard 
state in Sr2lr04 [l[ came as a surprise since this case 
departs from the established phenomenology of Mott- 
insulating behavior in the canonical early 3d transition- 
metal oxides. There, the localized nature of the 3d valence 
electrons is responsible for the small bandwidth W, large 
Coulomb repulsion U, and suppression of charge fluctua- 
tions In particular, Sr2ir04 appears to violate the 
U>W Mott criterion, which for the very delocalized 5d Ir 
electrons is not fulfilled. It was proposed that the strong 
spin-orbit (SO) interaction in 5d systems (Cso —485 meV 
for Ir [4]) might lead to instability against weak electron- 
electron correlation effects, and to the subsequent emer- 
gence of a many-body insulating ground state []]]. How- 
ever, the strong-SO limit J e ff — 1/2 ground-state scenario 
has recently been put into question la, and theoretical 
@ and time-resolved optical studies pj suggest that the 
insulating state of Sr2ir04 might be closer to a Slater 
than a Mott type: a band-like insulating state induced 
by the onset of antiferromagnetic (AF) ordering and con- 
sequent band folding at T/v — 240 K (Slater), as opposed 
to be driven by electron-correlations with an insulating 
gap already open at temperatures well above TV (Mott). 

Despite intense experimental and theoretical effort, the 
nature of the insulating state in the 5d iridates remains 
highly controversial. This is reminiscent of the situation 
in 3d oxides, for which the Mott versus band-insulator 
debate has lasted over four decades @, lil-fioj]. For in- 
stance, in the case of the prototypical AF insulator NiO, 
this debate was conclusively resolved only after the corre- 
lated nature of the insulating state was established based 
on: (i) the magnitude of the gap as measured by direct 
and inverse photoelectron spectroscopy (PES/IPES) fllj . 
much larger than expected from density functional theory 
(DFT) 8]; (ii) its persistence well above the Neel tem- 



perature T/v [12j; and (iii) the detailed comparison be- 
tween DMFT results [13j and momentum-resolved elec- 
tronic structure as measured by angle-resolved photoelec- 
tron spectroscopy (ARPES) fia [l4|. 

To address the nature of the insulating state in iri- 
dates, including the role of many-body electron correla- 
tions for their extended 5d orbitals and the delicate inter- 
play between W, U and SO energy scales, a particularly 
interesting system is the newly discovered AF insulator 
Na2lr03 [15|. Starting from a J e // = l/2 model in anal- 
ogy with Sr2lr04, this system was predicted to exhibit 
quantum spin Hall behavior, and was considered a po- 
tential candidate for a topologically insulating state [l6| ■ 
Further theoretical 17, 18 1 and experimental [l9[ work 
emphasized the relevance of structural distortions, which 
lower the local symmetry at the Ir site from octahedral 
(Oh) to trigonal (D^h)- Together with the structure com- 
prised of edge-sharing IrOe octahedra, this leads to an 
effective bandwidth for the Ir 5d-<2g manifold of ~ 1 eV. 
This potentially puts Na2lr03 closer than other iridates 
to the U ~ W Mott criterion borderline - and thus to a 
Mott insulating phase 17j. Most importantly, its lower 



Tv~15K provides the opportunity of studying the elec- 
tronic structure well above the long-range AF ordering 
temperature and - with the aid of novel DFT calcula- 
tions - establishing the nature of its insulating behavior. 

In this Letter we present a study of the low-energy elec- 
tronic structure of Na2lr03 by ARPES, angle-integrated 
PES with in-situ potassium doping, optics, and DFT 
calculations in the local-density approximation (LDA). 
The narrow bandwidth of the Ir-5d manifold observed in 
ARPES highlights the importance of SO coupling and 
structural distortions in driving the system towards a 
Mott transition. In addition, at variance with a Slater- 
type description, the gap is already open at 300 K and 
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FIG. 1. (color online), (a) Angle-integrated O and Ir valence- 
band photoemission spectrum of Na2lr03; the grey portion is 
shown in detail in (b-d). (b) ARPES EDCs for the Ir 5d-t 2g 
bands from along M—Y—K. (c) Negative second derivative of 
ARPES map highlighting the experimental dispersion; super- 
imposed (black diamonds) are the fit analysis results from (d). 
(d) Model fit of the F-point EDC with 4 Gaussian peaks for 
the Ir valence bands (VBs), and a Shirley background [20J: 4 
peaks are necessary to fit the dataset over the full momentum 
range [matching the number of DOS features from DFT in 
Fig.[4je,f)], with Gaussian lineshapes yielding better agree- 
ment than Lorentzians. Fhe average over M — Y — K gives 
(E vb ,YvbAEvb) in eV, for peak 1 to 4: (0.50,0.30,0.06); 
(0.94,0.39,0.08); (1.39,0.44,0.09); (1.89,0.43,0.07). 



does not show significant temperature dependence even 
across T/v — 15 K. From the potassium-induced chemi- 
cal potential shift and complementary optical conductiv- 
ity measurements, we estimate the insulating gap to be 



A, 



340meV. While LDA+SO already returns a de- 



pletion of density of states (DOS) at Ep, this only corre- 
sponds to a 'zero-gap'. The observed 340 meV gap value 
can be reproduced only in LDA+SO+U calculations, i.e. 
with the inclusion of both SO and U (with U ~ 3 eV) , es- 
tablishing Na2lr03 as a novel relativistic Mott insulator. 

The 130 K angle-integrated PES spectrum [2l[ in 
Fig.QTa) shows two broad spectral features belonging to 
the Ir 5d-t2g bands (0-3 eV binding energy), and to the 
O 2p manifold (beyond 3eV). The insulating character is 
evidenced by the lack of spectral weight at the chemical 
potential, which appears to be pinned to the top of the 
valence band (no temperature dependence is observed 



in the 130-250 K range [2ll |22|), Energy distribution 
curves (EDCs) measured by ARPES [H| along M-Y-K 
for the Ir 5d-t2 g bands are shown in Fig.QJb). The de- 
tected features are only weakly dispersing in energy, with 
the most obvious momentum-dependence being limited 
to their relative intensity. The electronic dispersion can 
be estimated from the negative second derivative map in 
Fig.UJc), calculated as — <9 2 /(k, E)/dE 2 , and more quan- 
titatively from the fit of EDCs as in Fig.QJd) (see cap- 
tion for details). The direct comparison of fit (black di- 
amonds) and second derivative results in Fig. [He) yields 
a good overall agreement in the dispersion of the 4 fea- 
tures (small deviations stem from the peaks' relative in- 
tensity variation, which is differently captured by the two 
methods). We find that the Ir 5d-t2 g valence band (VB) 
dispersions do not exceed AEvb ~ 100 meV in band- 
width - at variance with the generally expected larger 
hopping amplitude for bd-t2 g states. Another remarkable 
aspect of the Ir t^g bands is their linewidth, with values 
Yvb = V2ctvb ~ 300 — 450 meV. A possible origin might 
be many-body electron correlation effects as discussed for 
NiO [13j , and strong electron-phonon coupling leading to 
polaronic behavior in the spectral function [23j. 

The results in Fig.[T] already provide one very impor- 
tant clue: the gap is open well above Tjv — 15 K, which di- 
rectly excludes a Slater-type, magnetic-order-driven na- 
ture for the insulating state. As for the size of the gap, 
this cannot be readily identified by ARPES since photoe- 
mission can locate the valence band, as the first electron- 
removal state, but not the conduction band which be- 
longs to the electron-addition part of the spectral func- 
tion [2oj | . Alternatively, one can measure the gap in an 
optical experiment; however, one needs to discriminate 
between in-gap states of bosonic character (e.g., phonons, 
magnons, excitons) and those particle-hole excitations 
which instead determine the real charge gap. This com- 
plication can often hinder the precise identification of the 
gap edge [24[. Such complexity underlies the past con- 
troversy on NiO: while the 0.3 eV gap obtained by DFT 
Q was deemed consistent with optical experiments [25j |. 
the combination of direct and inverse PES revealed the 
actual gap value to be 4.3 eV [llj . This is well beyond the 
DFT estimate and establishes NiO as a correlated insu- 
lator. Here, for the most conclusive determination of the 
insulating gap magnitude, we use angle-integrated PES 
with in-situ doping by potassium deposition, as well as 
optics. A quantitative agreement between the two probes 
would provide validation against possible artefacts (26l |. 

To estimate the energy of the first electron-addition 
states and the DOS gap between valence and conduc- 
tion bands, we start by doping carriers (i.e., electrons) 
across the gap by in situ potassium deposition on the 
cleaved surfaces, and then follow the shift in chemical 
potential by angle-integrated PES. The results are sum- 
marized in Fig. [5] for K-evaporation performed at 130 K 
on two different freshly cleaved surfaces j22|. The most 
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FIG. 2. (color online), (a) Background-subtracted angle- 
integrated EDCs for selected values of K-exposure [see colored 
markers in (c)]; the chemical potential shift A/t is revealed 
by the motion of the high binding-energy trailing edge, (b) 
Same as in (a), but shifted by the corresponding Afi. (c) Afi 
vs. K-deposition time for two different samples; in the inset, 
K-induced low-energy spectral weight ASW as a function of 
A/t, for Sample 2 only. In panel (c) data and error bars are 
estimated from the comprehensive analysis of both O and Ir 
trailing edges and peak positions [2^]; also note that blue and 
red curves in (c) and its inset are both a guide-to-the-eye. 



evident effect is the shift towards higher binding energy 
of both Ir and O valence bands, as shown in Fig.[2Ja) 
for the Ir hd-ti g manifold [22j]. This arises from the 
(equal and opposite) shift of the chemical potential Afi 
when electrons donated by potassium are doped into the 
system; after an initial rapid increase, A/j, saturates at 
^340 me V [Fig.|2fc)]. When the K-deposited spectra are 
shifted in energy by the corresponding A/j, so that their 
high binding-energy trailing edges match the one of the 
fresh surface [Fig.^b)], one can observe the emergence 
of additional spectral weight (SW) in the region close 
to and above Ep. This low-energy K-induced spectral 
weight, ASW, can be computed as: 



ASW 



dk 



dE[I(k,E,x K +)-I(k,E,0)], (1) 



leV 



where I(k,oS) is the PES intensity, represents the 

K-induced surface doping, and Ep is the Fermi energy of 
the K-doped surface, which moves progressively beyond 
the undoped-surface Ep. The evolution of ASW plotted 
versus A/x in the inset of Fig.[5Jc) evidences an approxi- 
mately linear SW increase up to A/i ~ 200 — 250 meV, 



followed by a steeper rise once the saturation value 
A/i ~ 340 meV is being approached. This behavior can be 
understood as due to the initial filling in of in-gap defect 
states - either pre-existing or induced by K deposition - 
which makes the jump of the chemical potential not as 
sudden as for a clean insulating DOS. Only when elec- 
tronic states belonging to the Ir 5d-t2 g conduction band 
are reached one observes the saturation of A/j, and the 
more pronounced increase in ASW. This combined evo- 
lution of chemical potential shift and spectral weight in- 
crease points to a DOS insulating gap Ap£s~340meV. 

Turning now to the optical conductivity data (27j . in 
Fig. [3] we observe an insulating behavior with an absorp- 
tion edge starting at 300-400 meV at 300 K, with negligi- 
ble temperature dependence down to 8K and thus also 
across Tn — 15 K (see inset). We can fit the results us- 
ing a joint DOS with Gaussian peaks for the conduction 
band (CB) and each of the 4 VBs (see caption of Fig.QJ: 



4 r 

J(E) oc / dE'Ai G C b(E' + E) GvBi(E') 



(2) 



Here the prefactors A; represent the optical transition 
strengths, with the band index i running over the 4 VB 
features extracted from the ARPES data in Fig. [I] and 
arc left free. J(E) provides an excellent fit to the optical 
data in Fig. [31 and a least-squares analysis returns Eqb — 
680 meV for the location of the conduction band above 
Ep, with a width Tqb = V^ctcb — 160 meV (a value in 
agreement with DFT, as shown later). The consistency 
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FIG. 3. (Color online). Optical conductivity data (red line), 
together with the simulated Ir 5d-t2 g joint particle-hole DOS 
(black line), and its individual components from the simul- 
taneous fit of ARPES and optical data [colors and labels are 
consistent with the valence band features in Fig.[JJc), which 
represent the initial states of the lowest-energy optical tran- 
sitions]. Inset: temperature dependence of the gap edge. 
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FIG. 4. (color online), (a-c) Ir 5d-t2 g band structure (fc z = 0), and (d-f) corresponding DOS, obtained with LDA, LDA+SO, 
and LDA+SO+U (U — 3eV, J// = 0.6eV 30]). (g) k z dispersion, for the last occupied and first unoccupied Ir 5d-t2 9 bands from 
LDA+SO, as indicated by the filled region between the fc z =0 and n/c extreme lines. While the filled areas overlap, resulting 
in a vanishing indirect gap (i.e. Afc^O), the direct gap (i.e. Afc = 0) between valence and conduction bands is finite for all k z 
and fc|| . (h) LDA+SO direct gap distribution along M — V — K, for different fc z in the 3-dimensional Brillouin zone. 



of the combined PES-optical conductivity analysis is con- 
firmed by the optical gap obtained from the onset of the 
simulated conduction band, Aopt = E^J l g et . Following 
Ref.l23l. the latter is estimated as E^ et =E cb — Zocbi 
leading to Aqpt — 340 meV. This matches Ap^s from 
the K- induced Afi saturation in PES, providing a defini- 
tive estimate for the insulating gap, A gap ~340meV. 

With a gap much smaller than in typical Mott insu- 
lators, discriminating between correlated and band-like 
insulating behavior in NaalrOs requires a detailed com- 
parative DFT analysis [28j. Unlike the case of Sr 2 Ir04, 
the t 2g degeneracy and bandwidth in Na 2 Ir0 3 are af- 
fected by structural distortions and the presence of Na 
in the Ir plane. This is revealed by calculations we 
have performed for a distortion- and Na-free hypotheti- 
cal Ir02 parent compound: while in &O2 the individual 
^-dispersions are as wide as ~ 1.6 eV, the Na-induced 
band folding in Na2lrC>3 opens large band gaps, lead- 
ing to much narrower ^g-subbands (~ 100 meV). Even 
though this accounts well for the narrow bandwidth ob- 
served in ARPES, the material is still metallic in LDA 
with a high DOS at the Fermi level [Fig.|U[a,d)], at vari- 
ance with the experimental findings. Remarkably, when 
SO is switched on in LDA+SO a clear gap opens up at 
Ep in the k z = band dispersion [Fig.@|b)], although 
not in the DOS where only a zero-gap can be observed 
[Fig. (He)] . A closer inspection of the LDA+SO dispersion 
for first occupied and unoccupied bands - versus both k\\ 
and k z in Fig.[4jg) - reveals that the lack of a DOS gap 
stems from the overlap of VB and CB at T and K points 
for different k z values. In other words, while the direct 
gap (Afc = 0) is non-zero and ranges from a minimum of 
54 meV to a maximum of 220 meV over the full Brillouin 
zone [Fig.[4jh)], the indirect DOS gap (Afc ^ 0) is van- 
ishing. This is still in contrast with the experimentally 
determined band-gap magnitude A gap ~ 340 meV [29( . 



The disconnection between insulating behavior and on- 
set of AF ordering, together with the quantitative dis- 
agreement between observed and calculated gap even in 
LDA+SO, reveal that Na2lr03 cannot be regarded as ei- 
ther a Slater or a band insulator. Also, given the narrow 
tig bandwidths (~ 100 meV), one might expect the sys- 
tem to be even more unstable against local correlations 
than anticipated. Indeed, a good overall agreement with 
the data is found in LDA+SO+U, for U = 3eV and J H = 
0.6 eV [30j |: this returns a gap value A^Qg V ~ 365 meV 
[Fig.@ff)] close to the experimental A gap , and a 2eV en- 
ergy range for the Ir 5d-t2 g manifold [Fig.HJc)] match- 
ing the spectral weight distribution in Fig.Q] (note that 
a doubling of bands is seen in LDA+SO+U due to the 
imposed AF ordering, but is of no relevance here pjjl). 
At a first glance, U = 3 cV might seem a large value for 
bd orbitals; however, in the solid, the effective reduction 
of the atomic value of U strongly depends on the polar- 
izability of the surrounding medium, which is the result 
of many factors, in primis the anion-cation bond length 
[3~0 ] . In this perspective, the value we found is not unrea- 
sonable, and is also consistent with the existence of local 
moments above Tjy revealed by the Curie- Weiss magnetic 
susceptibility behavior with — 120K 15j . 



Our findings point to a Mott-like insulating state 
driven by the delicate interplay between W, U, and SO 
energy scales, in which co-participating structural distor- 
tions also play a crucial role. This establishes Na2lr03, 
and possibly other members of the iridate family, as 
a novel type of correlated insulator in which Coulomb 
(many-body) and relativistic (spin-orbit) effects cannot 
be decoupled, but must be treated on an equal footing. 
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